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ABSTRACT 


We have measured the reaction of O + Hjj' forming OH^ and H 2 O’''. This is 
one of the key gas-phase astrochemical processes initiating the formation of water 
molecules in dense molecular clouds. For this work, we have used a novel merged 
fast-beams apparatus which overlaps a beam of onto a beam of ground-term 
neutral O. Here, we present cross section data for forming and H 2 O’'' at rela¬ 
tive energies from 3.5 meV to ~ 15.5 and 0.13 eV, respectively. Measurements 
were performed for statistically populated O(^Pj) in the ground term reacting 
with hot HjI' (with an internal temperature of ~ 2500 — 3000 K). From these 
data, we have derived rate coefficients for translational temperatures from ^ 25 
K to ~ 10^ and 10^ K, respectively. Using state-of-the-art theoretical methods 
as a guide, we have converted these results to a thermal rate coefficient for form¬ 
ing either OH+ or H 20 ’'', thereby accounting for the temperature dependence of 
the O hne-structure levels. Our results are in good agreement with two inde¬ 
pendent flowing afterglow measurements at a temperature of ~ 300 K, and with 
a corresponding level of internal excitation. This good agreement strongly 
suggests that the internal excitation of the does not play a signihcant role 
in this reaction. The Langevin rate coefficient is in reasonable agreement with 
the experimental results at 10 K but a factor of ~ 2 larger at 300 K. The two 
published classical trajectory studies using quantum mechanical potential energy 
surfaces he a factor of ~ 1.5 above our experimental results over this 10 — 300 K 
range. 

Subject headings: astrobiology - astrochemistry - ISM; molecules - methods: 
laboratory - molecular data - molecular process 
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Introduction 


The genesis of life is believed to depend, in part, on the presence of water. Hence, 
nnderstanding interstellar formation of H 2 O is an important snbject for astrochemistry and 


astrobiology (Klippenstein et al. 2010). Gas-phase formation of water in both diffuse and 


dense molecular clouds is predicted to involve reactions of neutral O with H 3 (Smith & 


Spanel 1995) via 


0{^P) + B.tiK) ^ OH+ + H 2 , 
^ H 2 O+ + H. 


( 1 ) 

( 2 ) 


Through a series of subsequent hydrogen abstraction reactions with the abundant H 2 , the 
ionic products of Reactions ([^ and (|^ go on to form HaO"*". Dissociative recombination of 
HsO’*' with electrons results in several possible neutral products, one of which is H 2 O. The 
exact branching ratio for this specihc outgoing channel was uncertain for a long time, but 


has recently been measured to be ~ 20% (Novotny et al. 2010). At this point, one of the 


largest remaining uncertainties in the kinetics of the gas-phase formation of water appears 
to be the lack of reliable rate coefficients for Reactions ([^ and (|^. 

These two reactions also play a role in our understanding of diffuse and dense molecular 
clouds. The chemistry of these clouds is driven in part by the cosmic ray ionization rate 
of H 2 (CRIR) (. The temperatures of these environments are too cold for neutral-neutral 
reactions to overcome typical activation barriers. Instead, gas-phase chemistry proceeds 
largely by ion-neutral reactions ( Herbst fc Klemperer||1973 Smith fc Spanel||1995 ), a process 
that is initiated by cosmic ray ionization. The CRIR in diffuse clouds can be inferred using 
the OH+ and H 20 ’'' abundances. These are affected by Reactions ([^ and ([^ as described 


by Hollenbach et al. (2012). In dense clouds ( is constrained using the measured Hg 
abundance, which is also affected by Reactions ([^ and (|^ as is discussed in Klippenstein 


et al. (2010). Hence, constraining the astrophysical conditions in molecular clouds requires 
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reliable data for these two reactions. 

Published theoretical calculations for these reactions have been carried out using the 


classical Langevin method (e.g., Milligan & McEwan 2000). Classical trajectory studies 
using quantum mechanical potential energy surfaces (PESs) have also been published by 


Bettens et ah (1999) and Klippenstein et ah (2010). These latter two semi-classical results 


show reasonable agreement between one another. No fully quantum mechanical calculations 
exist as such calculations for reactions involving four or more atoms are too complex for 


current capabilities (Althorpe & Clary 2003; Bowman et al. 2011) 


Our previous experimental work for the analogous reactions (O’Connor et al. 2015) 

cCP) + ntCK) ^ CH+ + H2, 

^ CH+ + H, 


( 3 ) 

( 4 ) 


suggests that both the classical and semi-classical methods are likely to overestimate the 
rate coefficient for Reaction ([^, underestimate it for Reaction ([^, and incorrectly predict 
the temperature dependence of each. Unfortunately, the published room-temperature, 
flowing afterglow results for Reactions ([^ and ([^ cannot resolve this issue due to their 
large error bars ( Milligan fc McEwan||200^ and single temperature. Our work here aims to 
improve our understanding of the O -|- Hg reaction and address some of these outstanding 
issues. 


The rest of this paper is organized as follows. In Section we describe briefly our 
experimental apparatus and method. We present and discuss our results in Section 
and Section respectively. The astrochemical implications are discussed in Section A 
summary is given in Section 
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2. Experimental Apparatus and Method 


A detailed description of the ion-neutral merged-beams apparatus used for the present 


results can be found in O’Connor et ah (2015). Here we provide a brief description of the 
experiment and method, highlighting mainly those aspects which are new or specific to the 
present study. 


2.1. General 


Using a Cs-ion sputter source followed by a Wien filter, we formed a beam at a 
kinetic energy of Eq k. 28 keV (~ 1.75 keV/amu). This generated a pure beam of 0“(^P) 
as the oxygen anion possesses only a single bound term ( Rienstra-Kiracofe et al.||2002| ). The 
target material used in the source was AI 2 O 3 . 

The neutral O beam was produced via photodetachment of the anion inside a floating 
cell. The kinetic energy of the resulting O beam was controlled by varying the floating cell 


voltage t/f. The cross section for photodetachment of O has been measured by Lee fc 


Smith 


). The electron affinity of oxygen is 1.461 eV ( 

Rienstra-Kiracofe et al. 

2002 


Using our 808-nm (1.53-eV) laser with 1.8 kW of power, we estimated that approximately 
2% of the anions were converted to ground term Oi^P). The photon energy and number 
density were both too low to detach into higher terms. 


Similar to the work of Scheer et ah (1998) on photodetachment of C , we expect to 


statistically populate all three levels in the ground term of the atomic oxygen. The J = 1 
and 0 fine-structure levels lie above the J = 2 ground level by 19.6 meV and 28.1 meV, 
respectively. The thermal population of the J levels can be calculated using the partition 
functions 


5'je 


(-Ej/k^T) 




uj = 


( 5 ) 
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Here gfj = 2 J + 1 is the statistical weight of level J, k-Q is the Boltzmann constant, and T is 
the temperature. These populations are presented in Figure where one can see that they 
become statistical at temperatures above ~ 1000 K. In Section |4.2| we explain how we use 
our results with statistically populated ground-level O to generate thermal rate coefficients. 

The molecular beamline begins with a duoplasmatron source from which we extracted 
a beam of cations and used a Wien hlter to charge-to-mass select for Hj]’. The beam energy 
of E^+ ^ 5.29 keV (~ 1.75 keV/amu) was chosen to match the velocity of the oxygen anion 
beam. 


The formation mechanism of in a duoplasmatron leads to substantial internal 


excitation. In our earlier work of C on Hg (O’Connor et ah 2015), we inferred an internal 
temperature of ~ 2500 K. However, we found good agreement between our thermal rate 


coefficient results at 1000 K with the mass-scaled results of Savic et ah (2005), who studied 
C on Since the work of Savic et ah (2005) used DjI" with an internal temperature of 


77 K, the good agreement between their results and those of O’Connor et ah (2015) implies 
that the internal excitation of the H^ does not signihcantly influence reactions of the type 
X -|- HjI" —)■ XH^ -|- H 3 _„ for n = 1 or 2. We expect that this will also be the case for O on 


H^ and will return to this issue in Section 4.3 


In order to improve the H^ beam quality in the interaction region, we have modihed 
the beamline just before the cylindrical deflector used to merge the cation beam onto the 
neutral beam. In specihc, we installed a set of XY steerers and adjusted the location of 
the one-dimensional (ID) electrostatic lens before the merger. To compensate for the lack 
of vertical focusing in the cylindrical deflector, we used one pair of vertical steerer plates 
for focusing. One of the horizontal steering plates was used to adjust the angle of the 
beam going into the ID lens and cylindrical deflector. All of the other steerer plates were 
grounded. Additionally, the ID lens was moved closer to the cylindrical deflector so that we 
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could focus the beam at the appropriate location inside the deflector in order to generate a 
parallel beam at the exit (i.e., at the beginning of the interaction region). The resulting 
current in the interaction region was typically ~ 225 nA, corresponding to a typical number 
density of 10^ cm 

In the interaction region, the kinetic energy of any product ions formed was essentially 
the sum of Eq plus the product of the kinetic energy per amu (~ 1.75 keV/amu) times 
the mass in amu transferred from the H^. For O at 28 keV, the kinetic energy of the 
product was 29.75 keV for forming OH+, Reaction Q, and 31.5 keV for H20’'', Reaction 
(|^. At the end of the interaction region, the beams continued into a chicane which directed 
the Hg into a Faraday cup where the current was recorded. The daughter products were 
directed through the chicane and continue into an electrostatic energy analyzer. This 
consisted of a series of three 90° cylindrical deflectors with voltages optimized to direct 
the desired product ions into a channel electron multiplier (CEM). This hnal analyzer also 
served to discriminate against the dominant charged particle background which was due to 
the ~ 28 keV O’*" formed by stripping of the O beam on the residual gas downstream of 
the cation beam merger. The transmittance from the interaction region to the CEM was 
measured to be Ta = 0.74 ± 0.02. Here and throughout all uncertainties are quoted at an 
estimated Icr statistical conhdence level. 


2.2. Neutral Detector 


The neutral beam current Jq is monitored by measuring the secondary negative particle 


emission from a target inside a neutral particle detector (Bruhns et al. 2010a) and is given 
by 


Io = 


'ND 

w. 


( 6 ) 






Here Jnd is the current measured on the neutral detector, 7 is the secondary negative 
particle emission coefficient, and T„ is the transmission into the detector. Typical values for 
Iq are ~ 23 nA, as measured in amperes, with a statistical-like uncertainty of 5%. These 
currents correspond to particle number densities of 10 ^ cm 

We determined 7 using collisional detachment of 0“ on helium introduced into the 
chicane by a leak valve. This converts a portion of the initial 0“ beam to O and 0+. The 
positive and negative currents, Iq- and I 0 +, respectively, were measured in a Faraday 
cup, called the upper cup, which is situated behind a hole in the outer plate of the middle 
cylindrical deflector of the final analyzer. Because of the design of the final analyzer we 
can measure either Iq- or I 0 + in the upper cup simultaneous with Jnd, but not all three 
together. The transmittance from the interaction region into the upper cup was measured 
to be Tu = 0.64 ± 0.04. 


The total particle flux is assumed to be conserved at any given helium pressure p. This 
gives 


Jo_(p = 0) Io{p = 0) Io+{p = 0) = Io-{p) + loij)) + ^+(p), (7) 


where collisional detachment on the residual gas in the system generates non-zero currents 
of O and O’*" at p = 0. All of the currents are taken as positive quantities. In terms of 
measured quantities, we can rewrite Equation Q as 

ip-io) Jnd(o) /sno) _ /nd(p) i'^oAp) , , 

n r„ r„ r„ ’ ' ' 

where stands for the currents measured in the upper cup. Rearranging this yields 


—AJnd(p) 7u 


( 9 ) 


where A/(p) = /(p) — /(O). Typically we measure A/q_ and A/nd at set of a pressures 
Pi- Due to the coarse control on our leak valve, we measure AJq+ and AJnd at a slightly 
different set of pressures pk- Here we were able to match the pressures Pi and pk to 
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better than 3%. This small difference introdnces an insignificant nncertainty into onr 7 
determination. So we can convert these pk resnlts to those for pi nsing 

a/s*(p.) = a/s*(pJ^^^s?M (10) 

which we can then snbstitnte into Eqnation (|^. 

To determine 7 , we measnred /q_ and Jnd simnltaneonsly over intervals of ~ 250 s. 
First, no gas was introdnced in the chicane (p = 0). Next, we measnred with a helinm 
pressnre of Pi. Lastly, we measnred again withont helinm. For the cnrrents measnred 
withont gas, we nsed the average of the before and after measnrements. This allows ns 
to take into acconnt the finctnations of the 0“ beam dnring the measnrement. For each 
cnrrent measnred at high pressnre, we took the average of the cnrrent over the 250 s 
interval. We repeated this pattern to measnre Jnd and Iq+ at pk- Over the several months 
of onr O + HjI' measnrement campaign, we monitored 7 periodically and fonnd 7 = 2.6 ± 0.3. 
The nncertainty in 7 is treated as a systematic error. 


2.3. Beam Overlap and Relative Energies 


The overlap integral of the two beams has been determined nsing a combination of 
beam profile measnrements in the interaction region and geometric modeling, as described 
in 


Brnhns et al.| ( 

2010 b 

) and 

O’Connor et al. 

(2015 


between the two beams of 0.38 mrad was determined from the measnred beam profiles. 
This represents a factor-of-two rednction compared to onr C + H;]' work and is a resnlt of 
the modifications of the cation beamline described in Section 12.11 


The average relative energy (Ej-), along with the corresponding energy spread, was 
determined from the beam profiles measnrements and a Monte-Carlo simnlation of the 
beams trajectories described in O’Connor et ah (2015). Fignre shows the resnlts 
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values for {E.^) and the corresponding uncertainty. Based on these simulations, the 
lowest average relative energy achieved here is ~ 3.5 meV, corresponding to an effective 
translational temperature of ~ 27 K (as derived from a Maxwell-Boltzmann £t of the 
velocity distribution). We also find a Gaussian distribution for the interaction angle with a 


mean value of {6) = 0.71 ± 0.34 mrad. These improvements over O’Connor et al. (2015) are 
a direct result of the reduced bulk misalignment between the two beams. 


2.4. Merged-beams Rate Coefficient 


The measured merged-beams rate coefficient is the product of the cross section a and 
the relative velocity convolved with the energy spread of the experiment and is given by 

S e^VnVi 1 


(crUr) = 


( 11 ) 


IJ, L{n{z)). 

Here S is the signal count rate; Tg is the geometrical transmittance of the grid in front 
of the GEM; r] is the efficiency of the GEM; e is the elementary charge; and Uj are the 
laboratory velocities of the neutral and molecular ion beams, respectively; and h are the 
neutral and ion beam currents, respectively; L is the interaction region length; and {G(^)) 
is the overlap integral of the two beams in the interaction region. Typical values of all 
these parameters and their uncertainties are listed in Tables and Details about the 


associated uncertainties as well as the data acquisition procedure can be found in O’Gonnor 


et al. (2015). 


3. Experimental Results 


Our results for the merged-beams rate coefficient are presented in Figure as a 
function of {E^) for Reactions ([^ and ([^. Data were collected for merged-beam rate 
coefficient values down to ~ 1 x 10“^° cm^ s“^, below which the decreasing signal-to-noise 
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ratio made the required data acquisition times prohibitively long. From these, and 
using our calculated experimental energy spread, we can extract the cross section for 
each reaction for statistically populated ground-term O. We then generate translational 
temperature rate coefficients for each reaction by convolving the extracted cross section 
with a Maxwell-Boltzmann distribution describing the reaction center-of-mass velocity 
distribution. However, the internal energies of the reactants are not in thermal equilibrium. 


In Section |4.2| we describe how we can convert our translational temperature results into 
a thermal rate coefficient relevant for astrochemistry where the internal energies of the 
reactants are typically in thermal equilibrium. 


3.1. Cross Sections 


The cross section can be extracted from our data using the fitting function 

ao + aii2E^/'^ 


^ T. - 


( 12 ) 


Here x refers to either Reaction Q or Q, (Ja, is in units of cm^, and E is in eV. This 
function includes a term with an behavior at low energies that was chosen to match 

the predicted dependence in the thermal rate coefficient at low temperatures due 


to the charge-quadrupole interaction (Klippenstein et ah 2010). The other powers of E 


have been arbitrarily chosen to match the higher energy dependence in the measured 
merged-beams rate coefficients. 


To fit our measured data, we multiplied Equation (12) by and convolved the product 


with the experimental relative velocity distribution. The best fits to the data are shown in 
Figure by the solid lines. Tables gives the best fit parameters of the cross section for 
each reaction. The accuracy of the fit for Reaction Q is better than 7%. For Reaction ([^, 
the accuracy is better than 15%. Figure shows the experimentally derived cross sections 
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for each reaction. The results have been extrapolated to kinetic energies below ~ 3.5 meV 


using the theory of Klippenstein et ah (2010) as a guide. 


3.2. Translational Temperature Rate Coefficients 


We have derived the translational temperature rate coefficient for each reaction 
using the product of the extracted cross section times the relative velocity, all then 
convolved with a Maxwell-Boltzmann energy distribution. We fit the resulting translational 
temperature rate coefficient using 

Oo + + diT 




(13) 


ri/6 + (),/2 Ti/2 + bjT + 

where x refers to either Reaction 0 and ([^, ax is given in units of cm^ s and T in units 
of K. Table 1^ gives the best fit parameters for each reaction. 

Figure shows the experimentally derived translational temperature rate coefficients 
and associated uncertainties for Reactions ([^ and (|^. The uncertainty for each reaction, 
derived by adding in quadrature the estimated systematic error and the fitting error, is 
~ 14% for Reaction ([^ and ranges from 15 — 27% for Reaction (|^, for the respective 
energy ranges shown in Figure]^ The minimum value achieved for {Et) corresponds to an 
effective translational temperature of ~ 27 K, as discussed earlier. The fit function has been 
chosen so that the extrapolation of the translational temperature rate coefficient below 


27 K goes to a T dependence as predicted by Klippenstein et ah (2010). The highest 


values of {Et) measured for Reactions ([^ and (|^ of ~ 15.5 and 0.13 eV, respectively, 
yield approximate high temperature limits of ~ 10^ and 10^ K for the derived translational 
temperature rate coefficients. 
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4. Discussion 

Figure shows the asymptotic energy limits of various O + Hjj' reaction pathways. 
Reaction Q is exoergic by ~ 0.66 eV and Reaction (|^, by ~ 1.70 eV 
McEwan]|2000[ ). The energies are given for all parent and daughter products in their ground 
symmetries. 


(Milligan & 


4.1. Merged-beams Rate Coefficients 


4 . 1 . 1 . 0 + H+^0H+ + H2 

Taking into account the electronic states and spin symmetries of the reactants and 
products (Bettens et ah 1999), we can re-write Reaction ([^ as 


0(='F) + H+CA') ^ OH+(3s-) + H2(^S+). 


OH'*' can also be formed via the endoergic reaction: 


0(3f)+ H+(iA') ^ OH+(3 e-)+ H(2F)+ H(2F)-3.82eV. 


(14) 


(15) 


The energy dependence of our measured merged-beams rate coefficient for Reaction 0 
has a somewhat similar behavior to that measured for Reaction ([^ seen by O’Connor et al.| 


(2015). The merged-beams rate coefficient starts by increasing as the relative energy grows. 


This behavior could be due to the increasing number of ro-vibrational channels becoming 
energetically accessible in the reaction products or to the opening up of new electronic 
states in the intermediate reaction complex. Then, starting at about ~ 0.9 eV, the energy 
of the highest value of the merged-beams rate coefficient, there is a suggestion of a complex 


structure which was not seen in O’Connor et al. (2015). (All of the structures discussed here 


appears at lower energies than expected, by about 0.9 eV, as we explain in the following 
















- 14 


paragraph.) First the competing endoergic channel 


O + H+ ^ OH + H+ - 1.77eV 


(16) 


appears to open up, resulting in a reduction in the measured merged-beams rate coefficient. 
This is followed at ~ 3 eV by the opening up of another channel for forming OH'*', namely 


Reaction (15), causing a compensating increase in the signal. At higher energies, three 


additional competing channels open up: 


0 + H+ ^ O + H++ Ha - 4.34eV, 
^ 0++ H + Ha - 4.36eV, 
^ 0 + H + H+-6.16eV. 


(17) 

(18) 
(19) 


These are followed by a dramatic drop with increasing energy, which we attribute to the 
competing endoergic reaction pathways taking up all the flux of the reaction. 

The shift in these structures to energies lower than the known thresholds most likely 
comes from the molecular nature of the H^, with its range of possible ro-vibrational 


levels contributing to the process (as explained in more detail in O’Connor et ah 2015). 
This leads to a smearing out of the thresholds with the relative energy, blurring them 
together for the above-listed endoergic reactions. So we take the highest measured value 
of the merged-beams rate coefficient, at a relative energy of 0.89 eV, as the opening up of 


Reaction 16, which would be 1.77 eV for cold HA The lower-than-expected energy in the 


initial decrease of the merged-beams rate coefficient is probably due to the energy available 
from the internal excitation of H^ in our experiment. The ~ 0.87 eV energy difference lies 
within the ~ 0.5 — 1 eV range inferred from theoretical calculations (|Anicich fc Futr^ 


1984) and photodissociation measurements (X. Urbain, private communication). Using 


the Hg molecular partition function of Kylanpaa & Rantala (2011) and specihcally their 


Equation (8), this energy difference gives an internal temperature of ~ 3000 K. This is 
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slightly higher, but still consistent, with the internal excitation of ~ 2550 inferred in 


our C + Ht work which used the same Ht ion source (O’Connor et al. 2015). 


4.1.2. O + H+ ^ H 2 O+ + H 

Our merged-beams rate coefficient for Reaction ([^ decreases with increasing relative 
energy. Bettens et al. (1999) describe the formation of H20’'' as a two-step process. First 


a complex is formed where the oxygen atom extracts H+ from to form an OH+ ion 
adjacent to the remaining H 2 . Then the OH+ and H 2 undergo a relative rotation after 
which the OH'*' extracts a hydrogen atom from the H 2 . This entire process takes time and 
could explain the decrease seen in the rate coefficient. As the relative energy increases, 
the time available for the reaction also decreases. Clearly though, further theoretical and 
experimental studies are needed to better understand the process. 

At ~ 0.1 eV, we observe a strong decrease in the rate coefficient, similar to that of 
Reaction (|^ seen by O’Connor et al. (2015). The endoergic Reactions (15) - (19) are the 
probable explanation for this drop. The shift to a lower relative energy is also likely due to 
the internal excitation of H^. 


4.2. Converting Translational Temperature to Thermal Rate Coefficients 

In order to convert our translational temperature rate coefficients into a thermal rate 


coefficient, we follow the approach outlined in O’Connor et al. (2015). This method enables 


us to correct our results for statistically populated hne-structure levels in the atomic O 
into data relevant for thermally populated hne-structure levels. In this, we follow the 
approximations made by Klippenstein et al. (2010), in their calculations for the O -|- Hjj" 
collision system, namely that the reaction proceeds adiabatically and that surface crossings 
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and intersystem transitions are both nnimportant. 

Another assnmption that we make is that the internal excitation of the Hjj’ does not 
affect the reaction. Theoretically, this approach is snpported by the long-range character of 
the entrance PES being dictated by the polarization of the hne-strnctnre state of the oxygen 


in the held of a point charge (Gentry & Giese 1977). Experimentally, it is snpported by the 


good agreement between our G -|- results on internally excited Hjj' and the mass-scaled 
results of Savic et al. (2005) for G on D;^ with an internal temperature of 77 K. It is also 
supported by the good agreement that we hnd below between our thermal O -|- results 


and those of Fehsenfeld (1976) and Milligan & McEwan (2000) for Hg with an internal 
termperature of 300 K. 

Taken all together, these approximations enable us to derive temperature-dependent 
multiplicative scaling factors to convert our translational temperature results into thermal 


results. For this we use the theoretical study of Gentry & Giese (1977) for ground 


term atomic oxygen in the presence of a positive charge. Bettens & Gollins (2001) and 


Klippenstein et al. (2010) also based their calculations on Gentry & Giese (1977), which 


found that the nine states in the 0 (^P) manifold are split into three attractive and 
six repulsive surfaces at long-range separation of the reactants. The Mj = 0 and ±1 
components of the hve-fold degenerate ground state ^P 2 state correlate with the attractive 
PES, while the Mj = ±2 components of the ^P 2 and all components of ^Pi and 
correlate with the repulsive ^fl surface. The partition functions for the attractive and 
-reactive ^fl PES shown in Figure are given by 

( 20 ) 


non- 


Ms = yM 2 
5 

2 

Mn — Mo + Ml -|- —U2 

5 


( 21 ) 


with Mj dehned by Equation (|^. 


The theoretical approach outlined above provides information only for the total 
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reaction probability, i.e., the sum of Reactions Q and ([^. We are aware of only one 
published experimental study for the branching ratios for forming OH’*’ and H 2 O’'’, 


namely the results of Milligan & McEwan (2000) at 300 K. How those results scale with 


temperature is uncertain. Indeed, O’Connor et ah (2015) found a signihcant temperature 
dependence for the branching ratio in the analogous C + Hjj’ reaction. Lacking the necessary 
temperature-dependent information on the branching ratios, in order to convert our 
translational temperature results to a thermal rate coefficient, we sum our translational 
temperature rate coefficients for Reactions ([^ and (|^ and then multiply the result by 3u^. 
This factor of 3 takes into account the fact that in our measurement only one-third of the 
O hne-structure levels contribute to the reaction process. The resulting summed thermal 
rate coefficient is shown in Figure 


We can also use this theoretical approach to generate a temperature-dependent 
Langevin rate coefficient. The classical Langevin rate coefficient of 1.3 x 10“® cm^ s“^ 

, considers all of the O -|- H(|’ symmetries involved in the reaction 
process to be attractive. We can convert this value to a temperature-dependent thermal 
rate coefficient by multiplying it by the partition function us. Both the unmodihed and 
modihed Langevin rate coefficient are shown in Figure]^ 


(Milligan fc McEwan|2000) 


4.3. Thermal Rate Coefficients 

Our experimentally derived thermal rate coefficient for O -|- Hjj’ forming either OH’*’ 
or H 2 O’*’ is in good agreement with the unmodified Langevin rate coefficient at 10 K. 
At this temperature, the modified Langevin rate coefficient lies a factor of ~ 1.6 below 
our experimental results. This is not surprising since the charge-quadrupole interaction 
responsible for the growth of the rate coefficient at low temperatures is absent from the 
Langevin model. As the temperature increases, the unmodified rate coefficient becomes 
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increasingly discrepant with the experimental results, growing to a factor of ~ 2 times 
larger at 300 K. Conversely the modihed rate coefficient comes into reasonable agreement 
with our results, lying just outside the experimental error bars, but matching the general 
temperature dependence. 


The semi-classical results of Bettens et ah (1999) and those of Klippenstein et ah 


(2010) both shown in Figure]^ each display a trend similar to our experimentally derived 
thermal rate coefficient. However, both sets of values lie outside our experimental error 
bars. The Bettens & Collins results are a factor of ~ 1.7 larger at 10 K and a factor of ~ 


1.3 at 300 K, while the Klippenstein et ah (2010) results are a factor of ~ 1.4 and 1.7 larger 


at 10 and 300 K, respectively. The cause for the discrepancy is not clear. For the analogous 
summed thermal rate coefficient for C H^, reasonable agreement was found between our 


results and the calculations of Bettens & Collins (1998, 2001) 


Our derived thermal rate coefficient is in good agreement with the experimental results 


from Fehsenfeld (1976) at 300 K. In his flowing afterglow experiment, he was unable 


to distinguish the products of the reaction and thus gave only an overall thermal rate 


coefficient for O on Hg . Milligan & McEwan (2000) measured the rate coefficients of 
Reactions ([^ and (|^ at 295 K. We compare with their summed results and find that their 
error bar overlaps with ours. So even though their results are a factor of ~ 1.8 times larger 
than ours, the agreement seems reasonable enough, given the challenge of monitoring the 
atomic oxygen density in the flowing afterglow method. 


Comparing all of the experimental thermal results for O on at ~ 300 K, we find good 
agreement between our work with hot (~ 3000 K) and the two published measurements 
on cold 300 K). Similarly good agreement was found at a temperature of ~ 1000 K 


between our C -|- Hg work compared to the mass-scaled results of Savic et ah (2005) for 


C -I- D^, which had an internal temperature of 77 K. Taken together, these experimental 
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results strongly suggest that for reactions of the type X + the internal excitation of the 
does not play a signihcant role in determining the total rate coefficient for reacting, 
summed over all possible outgoing channels. We hypothesize that this behavior is due, in 
part, to the reaction being driven primarily by the long range portion of the PESs for the 
reaction systems, where the internal excitation of the molecule is not yet felt. Furthermore, 
has only a single bound electronic symmetry and so the entire population has the same 
spin symmetry. Thus any spin conservation in the reaction affects the entire population 
the same way, as opposed to a molecule with a population divided between two or more 
different spin multiplicities. Clearly, though, further experimental and theoretical work is 
needed to more solidly understand the physics of the reaction. 


5. Astrochemical Implications 


We have used the gas-phase astrochemical code Nahoon (Wakelam et ah 2012) and the 


2014 version of the Kinetic Database for Astrochemistry (KIDA; Wakelam et al. 2015) in 
order to investigate some of the astrophysical implications of our Endings, particularly for 
the gas-phase chemistry of dense molecular clouds. Though a complete study is beyond 
the scope of this paper (for example, we leave for future study the eEects of surface or 
ice chemistry as well as any depletion efiects), these preliminary models do give some 


insight into the impact of our results. Following Wakelam et al. (2015), the specific 
cloud parameters used were a hydrogen nuclei density of uh = 2 x 10^ cm“^, a CRIR of 
C = 10“^^ s“^, and a visual extinction of Ay = 30. We also adopted the initial abundances 
given in their paper. 


For these simulations we have taken our experimentally derived thermal rate coefficient 
for O -I- HjI' for forming either OH+ or H20^ and implemented it into Nahoon/KIDA. 
Since the available fitting functions within Nahoon/KIDA were unable to accurately fit our 
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thermal rate coefiiient, we have ht it using the recommended function of Novotny et al, 


(2013), namely 


k{T) = A 


300 K 
T 


+ r-=>/‘"^Ciexp(-ri/r). 


( 22 ) 


2=1 


The resulting £t parameters are given in Table The deviations of the fitted k from the 
data are less than ~ 0.5% over the reported 1 — 10^ K temperature range. 


We have implemented our thermal rate coefficient into Nahoon/KIDA, modifying 
Nahoon to handle the above equation. Since we do not know the branching ratios for the 
formation of OH'*' and H 20 ’*', we have run simulations hrst assuming a branching ratio of 
100% for forming OH”*" and later assuming 100% for forming H 20 ’''. This is justified as 
both ions are predicted to be highly transitory, rapidly reacting with the abundant H 2 in 
the cloud to form HaO’*' via hydrogen abstraction, 

OH+ -% H 2 O+ H 3 O+. (23) 

As we expected, we found no signihcant difference in our model results for either branching 
ratio assumption. 


5.1. Formation of Water 

To investigate the impact of our results on the formation of water, we have calculated 
the HaO’*' and H 2 O abundances for a dense molecular cloud at a temperature of 10 K. Over 
most of the lifetime of a cloud, HaO"*" is the dominant gas-phase precursor of water, forming 
H 2 O via dissociative recombination 

e" + H 3 O+ ^ H 2 O + H. (24) 


Figure shows the ratios of the predicted abundances using our results relative to those 
from the unmodihed Nahoon/KIDA, which uses the theoretical rate coefficient results of 
Klippenstein et al. ( 2010| . 
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At a cloud age of 10^ years, the abundance of both HsO’*' and H 2 O are reduced by 
~ 15%. These changes can be traced back to our reduced rate coefficient for forming OH’*' 
or H 20 ''', which go on to HsO"*" via hydrogen abstraction and then water via dissociative 
recombination. Naively, one might expect the HsO’*' and H 2 O abundances to be reduced by 


the ratio of our rate coefficient relative to that of Klippenstein et ah (2010), which is ~ 0.72 


at 10 K. However, our lower rate coefficient results in a reduced destruction rate for the 
and thus a higher predicted abundance, giving an abundance ratio at 10 ^ years of ~ 1.18 
compared to the unmodihed Nahoon/KIDA results. Multiplying these two factors together 
explains the ~ 15% reduction seen at short times. 

At longer times, the observed reductions in HsO^ and H 2 O become much more 
signihcant. Identifying the cause for this, though, is complicated, due to an increase 
in the complexity of the relevant chemical network. As an example, we discuss here 
the intermediate times between ~ 10^'® and 10^'® years. During this time, along with 


reaction (24), two additional important channels for forming water open up, namely 


C3 + HqO^ — y H2O + 


(25) 


and 

HCN + H 3 O+ ^ H 2 O + HCNH+. (26) 

Indeed, these two channels are predicted to dominate over dissociative recombination for 
part of this epoch. At longer times, dissociative recombination again comes to dominate. 

New channels for forming H 30 ’'' also open up between ~ 10"^'® and 10^ ® years, including 

H 2 O + H+ ^ H 2 + H 3 O+ (27) 


and 


H2O + HCO+ ^ CO + H3O+. 


( 28 ) 
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However, these reactions lead to destruction of H 2 O as well, since only a fraction of the 
resulting HsO"*" gets converted back into water. We hnd that with the modihed chemistry, 
the larger decrease seen in the relative H 2 O abundance (compared to that for HsO^) can 
be traced back, in part, to the increased Hg abundance during this epoch, leading to an 
enhancement in HCO'*' formation via 


CO + H+ ^ HCO+ + H 2 . 


(29) 


The resulting HCO’*' then leads to an increased reduction in the water abundance via 


reaction (28) for the modified chemistry. 


5.2. Cosmic Ray Ionization Rate 

In dense molecular clouds, the CRIR ( can be constrained through observations of 


the Hg column density (e.g., McCall et ah 1999 Klippenstein et ah 2010 Oka 2013). 
The H+ is formed in a two-step process, beginning with the ionization of H 2 by cosmic 
rays; other ionization processes of H 2 are insignihcant. At typical dense cloud densities of 
nn ~ 10^ cm“^, the resulting rapidly reacts with a neutral H 2 molecule, in less than a 
day dOka 2013), meaning that generation is completely determined by the CRIR, with 
a formation rate versus time t of 

'dn. 


''Ht 


dt 


C^H2 5 


(30) 


form 


where nHj is the number density of molecular hydrogen. The dominant destruction 
mechanisms of in dense clouds are binary chemical reactions with other species, 
primarily neutral since the electron abundance is too low for dissociative recombination to 
play a significant role. Thus we can write the destruction rate as 


dnH+ 

df 


dest 




Tiin. 


(31) 
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where ki is the thermal rate coefficient for the reaction of species i with and rii is the 
number density of species i. 


In the quasi-equilibrium of dense molecular clouds, we can equate Equations (30) and 

Cn^2 = ^kinin^+. (32) 


(31) giving 


As described in McCall et ah (1999), Klippenstein et ah (2010), and Oka (2013), the 


properties of dense clouds enable us to make the approximation 




(33) 


where N-^+ is the observed column density of Hj)" and L the absorption path length through 


the cloud. Combining Equations (|3^ and (pl3|) and rearranging gives 

J2ki 


CL = 




UHo 


iT^i 


(34) 


The quantity C,L can thus be constrained by a combination of the observed column 
density and riHa number density and astrochemical calculations of the individual terms 


contributing to the summation over species. Previous studies into Equation (32) have 


limited the number of species considered. For example, McCall et ah (1999) considered 


only reactions with CO, Klippenstein et ah (2010) also included O, and Oka (2013) added 


to these N 2 and electrons. Here we have used the modified Nahoon/KIDA to calculate 
^ • kiTii for the most significant species as well as for all species. For comparison, we have 


also calculated the unmodified Nahoon/KIDA results for all species. Figures 10 and 11 
respectively, show our results for cloud ages of 10^ and 10® years over the temperature range 
of 10 — 400 K. These times lie within the range of commonly inferred cloud ages. 


At 10® years and below ~ 160 K, the most important reactants with H;/ are CO, O, 
C, H 2 O, and the HCN/HNC isomers, in descending order of importance. The contributions 
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of each of these species are shown by the various colored solid lines in Fig. 10 and the 
total ^i'^i due to all species by the solid black line. Above ~ 250 K, the abundance of 
C decreases dramatically and reactions with C become unimportant. It is also interesting 
to note that the structure due to CO largely gets washed out by the contributions of O, 

C, and H 2 O. Additionally, the structures seen in the total fcjUj are due primarily to 
the HCN/HNC isomers between ~ 120 and 250 K and due to water above 250 K. The 
increase due to H 2 O is a result of the higher temperatures enabling neutral-neutral reactions 
leading to water to become important. Using the unmodihed Nahoon/KIDA yields a total 
summation shown by the black dashed line, which is up to ~ 15% larger than the modihed 
value. 


At 10® years and below ~ 240 K the most important reactants with are CO, O, 
and N 2 , in decreasing order of importance. The contributions of each of these species 


are shown by the various colored solid lines in Fig. 11 and the total due to all 

species by the solid black line. Above ~ 240 K, the abundance of O decreases dramatically 
and reactions with O become unimportant. Conversely, the abundance of H 2 O increases 
as neutral-neutral reactions leading to water become important. Using the unmodihed 
Nahoon/KIDA yields a total summation shown by the black dashed line, which is up to 
larger than the modihed value. 


These hndings indicate that constraints to the value (L depend, in part, on using a 
complete chemical model, knowing the age and temperature of the observed cloud, and 
using accurate reaction rate coefficients. The calculated abundances rii also depend on the 
adopted values for C, and riHa astrochemical model, implying that one will have to 

iterate the model in order to achieve convergence with the observations for the quantity 
C,L. Our experimentally derived thermal rate coefficient reported here helps to improve the 
reliability of this approach to determining (/L, but it is clear that there are many additional 




parameters which also need to be well constrained. 


6. Summary 

Using a novel merged fast-beams apparatus, we have measured the cross sections 
for O -|- Hjj' forming OH’*' and H20’''. Our measurements were performed for statistically 
populated O(^Pj) in the ground term and internally hot (~ 2500 — 3000 K). Using 
state-of-the-art theory as a guide to account for the temperature dependence of the O 
hue-structure J-levels, we have converted our results into a thermal rate coefficient for 
forming either OH'*' or H20^. The good agreement that we hud with the two published 
howing afterglow meaurements at a temperature of ~ 300 K (and a corresponding level 
of Hjj' internal excitation) strongly suggests that the internal excitation does not 
signihcantly affect the thermal rate coefficient for this reaction. The Langevin value is 
in good agreement with our results at 10 K but a factor of 2 higher at 300 K. The two 
published semiclassical results lie a factor of ~ 1.5 above our results over this temperature 
range. We have implemented our results into the astrochemical code Nahoon to explore 
some of the astrophysical implications of our results. For example, for dense clouds at 10 K, 
we hnd a reduction in the predicted water abundance by up to nearly 40% at certain times 
in the lifetime of the cloud. 
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0 12 3 4 

|C/f| (kV) 


Fig. 2.— Simulated average relative energy as a function of the floating cell voltage 
\Uf\. The vertical error bars show the full width at half maximum spread of the calculated 


energy distribution. 
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r-; 



Fig. 3.— Experimental rate coefficients {(tv-^) as a function of the average relative energy 
(Er) are shown for Reaction Q by the filled circles and for Reaction ([^ by the open squares. 
The error bars signify the Icr statistical uncertainties. Empirical fits to the experimental data 


using Equation (12) are shown by the solid lines. The vertical lines at the top of the graph 


show the threshold energies for the competing Reactions (16)-(19) of 1.77, 4.34, 4.36 and 
6.16 eV, respectively. 
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Relative energy (eV) 


Fig. 4.— The solid black lines show the experimentally derived cross section a as a function 
of the average relative energy {E^) for Reactions Q and (|^. The quadrature sum of the 
13% systematic uncertainty and the fitting accuracy is denoted by the shaded region. The 
extrapolation of the experimental results to lower kinetic energies is shown by the red dashed 
lines and the systematic uncertainty, shown by the surrounding shaded region, is taken as 
constant using that of the lowest measured energy. 
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Fig. 5.— Experimentally derived translational temperature rate coefficients vs. tempera¬ 
ture for Reactions Q and ([^ are shown by the black solid lines. The quadrature sum of 
the systematic uncertainty and the htting accuracy is denoted by the shaded region. The 
red dashed lines are extrapolations to lower temperatures with the systematic uncertainty 
(shaded region) taken as constant using that of the lowest measured energy. 
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Fig. 6.— Energy level diagram for various O + reaction pathways, given in eV, for the 
various systems in their ground state. Values for HaO^, H20’'' + H and OH^ + H 2 are given 


by Milligan & McEwan (2000); the others are derived from Huber & Herzberg (1979), Rohse 


et al. 

(1994) 

Liu et al. 

(2009 

) and 

Kramida et al. 

(2014 
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Temperature (K) 


Fig. 8.— Our experimentally derived thermal rate coefficient vs. temperature for the sum of 
Reactions Q and ([^ is shown by the black solid line. The dashed lines are the Langevin rate 
coefficient unmodihed (thin) and modihed (thick). The dotted-dashed lines present the total 


theoretical calculations of Bettens et ah (1999). The Klippenstein et al. (2010) calculations 
are shown by the dotted line. The inverted open triangle is the 300 K experimental result 


of Fehsenfeld (1976) shifted to 330 K for clarity, and the hlled triangle the total 295 K 


experimental result of Milligan & McEwan (2000). 























36 



Fig. 9.— Ratio of the modified to unmodified predicted abundances for HsO’*' (dashed line) 
and H 2 O (solid line). The ratio gives the modihed Nahoon/KIDA results using our exper¬ 
imentally derived total rate coefficient for the O -|- Hjj" reaction divided by the unmodihed 


Nahoon/KIDA results, which uses the rate coefficient from Klippenstein et ah (2010). 
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Fig. 10.— Predicted factor - kiUi at a cloud age of 10® years using the modified and 
unmodified Nahoon/KIDA. The solid colored lines correspond to the number of terms being 
included in the summation for the modihed factor, starting with dark blue for just CO, red 
for the addition of O, green for C, light blue for H 2 O, violet for HCN/HNC, and black for 
the addition of all other species. The predicted factor for all reactants using the unmodihed 
Nahoon/KIDA is shown by the black dashed curve. 




ox) 
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Fig. 11.— Same as Fig. 10 but for a cloud age of 10° years. Here the cumulative additions of 


the reactants i = {CO, O, N 2 , H 2 O, all other species} are shown by the dark blue, red, green, 
light blue, and black solid lines, respectively. 
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Table 1: Typical values of the statistical-like uncertainties for Equation 0 for a single 


data run. 


Source 

Symbol 

Value 

Uncertainty (%) 

Signal rate 

5 

1-15 Hz 

< 4 

0 velocity 

Vn 

5.8 X 10^ cm s“^ 

<1 

Hg velocity 

Vi 

5.8 X 10^ cm s“^ 

<1 

0 current 

In 

23 nA 

5 

Hjj' current 

li 

225 nA 

5 

Overlap factor 

{n{z)) 

2.7 cm“^ 

10 

Statistical-like uncertainty (single run) 



13 


Note. - The confidence level for each value is taken as equivalent to a Icr confidence level. 
The total statistical-like uncertainty is calculated by treating each individual uncertainty as 
a random sign error and adding all in quadrature. 
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Table 2: Same as Table 1 but for the systematic uncertainties for all data runs. 


Source 

Symbol 

Value 

Uncertainty (%) 

Analyzer transmission 

T 

0.74 

3 

Grid transmission 

T 

0.90 

1 

Neutral transmission 

T 

0.94 

2 

Neutral detector calibration 

7 

2.6 

12 

GEM efficiency 

7] 

0.99 

3 

Interaction length 

L 

121.5 cm 

2 

Total systematic uncertainty 



13 
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Table 
in eV, 


3: Fit parameters for the cross section of Reactions ([^ and Q in units of cm^ for E 


using Equation (12). 


Reaction 


( 1 ) 

( 2 ) 


Parameters 


Oq 0 - 1/2 ^2 ^4 

3.7314E-16 2.1237E-16 - 3.6414E-02 5.0532E-04 

9.8531E-17 - -4.0668E-02 -4.1891 517.80 
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Table 4: Fit parameters for the 
in units of cm^ s“^ for T in K, 


kinetic temperature rate coefficient for Reactions 0 and (|^ 


using Equation (13). 


Reaction 




Parameters 






Oq 

fll/2 

Oi 

^1/2 

hi 

^3/2 

( 

1 

) 

5.1142E-10 

2.6568E-11 

9.8503E-15 

1.6747E-02 

-9.9613E-5 

1.1006E-6 

( 

2 

) 

4.2253E-10 

- 

- 

3.4977E-03 

-1.4126E-04 

6.3584E-05 











Table 5. Fit parameters using Equation (22) for the thermal rate coefficient of O + H 


forming either or H 20 ’^. 


Parameter 

Value 

X 

y 

Units 

A 

7.39 - 

-10 

S —1 

s 

n 

1.46 

-1 

dimensionless 

Cl 

6.32 

-8 

K3/2 g-i 

C2 

-3.00 

-6 

K3/2 g-1 

C3 

-1.17 

-5 

K3/2 g-1 

C4 

5.76 

-4 

K3/2 g-1 

Ti 

7.47 

1 

K 

T2 

5.40 

2 

K 

Ts 

1.83 

3 

K 

T4 

1.90 

4 

K 


Note. — The value for each parameter is 


given by X X 10^. 








